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Motivated by the Kronecker product approximation technique, we have developed a very 
simple method to assess the inseparability of bipartite quantum systems, which is based 
on a realigned matrix constructed from the density matrix. For any separable state, 
the sum of the singular values of the matrix should be less than or equal to 1. This 
condition provides a very simple, computable necessary criterion for separability, and 
shows powerful ability to identify most bound entangled states discussed in the literature. 
As a byproduct of the criterion, we give an estimate for the degree of entanglement of 
the quantum state. 
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1. Introduction 

Quantum entangled states have recently emerged as basic resources in the rapidly expand- 
ing field of quantum information processing, with remarkable applications such as quantum 
teleportation, cryptography, dense coding and parallel computation ^ |2J E] ■ However, two 
fundamental questions have only been partially answered: how do we know if a given quantum 
state is entangled, and how entangled is it still after interacting with a noisy environment? 

From a practical point of view, the state of a composite quantum system is called unentan- 
gled or separable if it can be prepared in a " local" or "classical" way. A separable bipartite 
system can be expressed as an ensemble realization of pure product states \ip%}^ \4>i) b occur- 
ring with a certain probability pi : 

PAB =^,Pipt ® Pi \ (!) 

i 

where pf = \ipi) A (ipi\, pf = \<f>i) B (<^|, J2iPi = !> and IV'i),,!' I^*)b are normalized pure 
states of subsystems A and B, respectively 0]. If no convex linear combination exists for a 
given pab, then the state is called " entangled" . 

There have been considerable efforts in recent years to analyze the separability and 
quantitative character of quantum entanglement. For a pure state pab, it is separable iff 
Pab = PA <8> Pb where pa,b is the reduced density matrix defined as pa — TtbPab and 
Pb = TtaPab- However, for a generic mixed state pab, finding a decomposition as in Eq. 
or proving that it does not exist is a non-trivial task (see |Sj and references therein). The 
first important breakthrough was made by Peres who proposed that partial transposition 
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with respect to one subsystem of the density matrix for a separable state is positive, i.e., has 
non-negative eigenvalues [HJ ■ Now known as the PPT criterion, this was shown by Horodecki 
et al to be sufficient for bipartite systems of 2 x 2 and 2 x 3 jZj. In the same paper, a necessary 
and sufficient condition for separability was found by establishing a close connection between 
positive map theory and separability. Soon after, Wootters succeeded in computing the "en- 
tanglement of formation" jH) and thus obtained a separability criterion for 2x2 mixtures [HJ ■ 
The "reduction criterion" proposed independently in and gives another necessary cri- 
terion which is equivalent to the PPT criterion for 2 x n composite systems but is generally 
weaker than the PPT criterion. Pittenger et al gave also a sufficient criterion for separability 
connected with the Fourier representations of density matrices ^21- Later, Nielsen et al [T3] 
presented another necessary criterion called the majorization criterion: the decreasingly or- 
dered vector of the eigenvalues for pab is majorized by that of pa or pb alone for a separable 
state. A new method of constructing entanglement witnesses for detecting entanglement was 
given in |7| and [141 115| . There are also some necessary and sufficient criteria of separability 
for low rank cases of the density matrix, as shown in [161 117) . In addition, it was shown in 
P| and Q2| that a necessary and sufficient separability criterion is also equivalent to certain 
sets of equations. 

However, despite these advances, a practical computable criterion for generic bipartite 
systems is mainly limited to the PPT, reduction and majorization criteria, as well as a recent 
extension of the PPT criterion based on semidefinite programs [201 - The PPT criterion has 
been considered a strong one up to now, but in general it is still not sufficient for higher 
dimensions of greater than or equivalent to 3. Several counterexamples of bound entangled 
states with PPT properties were provided in |21|. and their entanglement does not seem to 
be "useful" for distillation (22]. These states are "weakly" inseparable and it is very hard to 
establish with certainty their inseparability jllj . 

In this paper we focus on an inseparability test for a generic bipartite quantum system of 
arbitrary dimensions. The basic procedure is the same as that of the cross norm criterion pre- 
sented in but is mathematically much more straightforward. Motivated by the Kronecker 
product approximation of the density matrix [241 125) we have derived a directly computa- 
tional method to recognize entangled states based on a realigned matrix constructed from the 
density matrix. An estimate is also given for depicting the "distance" from a quantum state 
to the maximally entangled state and to the maximally mixed state. Several typical examples 
of bound entangled states are shown to be recognized by this criterion in Section 3. A brief 
summary and discussion are given in the last section. 

2. A necessary separability criterion based on a matrix realignment method 

In this section we will present an inseparability criterion to recognize entangled states based 
on simple matrix analysis. Some of its characteristics are shown in several propositions. At 
the same time, we give a possible measure for the degree of entanglement as a by-product. The 
main tools used are a technique developed by Loan and Pitsianis |2H] for the Kronecker 
product approximation of a given matrix, and some results from matrix analysis (see Chapters 
3 and 4 of |2(i| for a more extensive background). We shall first review some of the notation 



Kai Chen and Ling-An Wu 3 



and results required in this paper. 

Definition: For each m x n matrix A — [a%j\, where a%j is the matrix entry of A, we define 
the vector vec(A) as 

vec(A) = [an, • • • , a m i, 0,12, • • • , a m 2, • • • , ai n , ■ • • , a mn ] T . 



Let Z be an m x m block matrix with block size n x n. We define a realigned matrix Z of 
size m 2 x n 2 that contains the same elements as Z but in different positions as 

vec(Zi. 1 ) T 



Z 



vec(Z mA ) T 
vec(Zi tm ) T 
vec(Z m ^ m ) T 

so that the singular value decomposition for Z is 

Z = UT,V^ = °-, 



(2) 



II- iV 



(3) 



where U — [uiu 2 ---u m 2] g C m xm and V = [viv 2 ---v n 2] e C n x " are unitary, S is a 
diagonal matrix with elements <j\ > (T2 > • ■ ■ ~> & q > and g = min(rn 2 ,n 2 ). In fact, 
the number of nonzero singular values Oi is the rank r of matrix Z 1 and <Ji are exactly 
the nonnegative square roots of the eigenvalues of ZZ^ or Z^ Z 26 . Based on the above 
constructions, Loan and Pitsianis obtained the following representation for Z 



Z = ® Yi), 



(4) 



i=i 



with vec(Xi) = y/alui and vec(Yi) = y/oiv* |24l 125) . 

For any given density matrix pab we can associate a realigned version pab according 
to the transformation of Eq. J2J. For example, a 2 x 2 bipartite density matrix p can be 
transformed as: 



/ P11 P12 

p2l P22 



Pai P32 
\ Pai PA2 



Pi3 pu \ 

P23 P24 



P33 P34 
P43 P44 / 



P = 



Pll 


p2l 


P12 


P22 


P31 


Pa 


P32 


P42 


Pl3 


P23 


Pu 


P24 



(5) 



\ P33 P43 P34 P44 / 

For any separable system pab, we arrive at the following main theorem 



Theorem: If an mx n bipartite density matrix pab is separable, then for the m 2 x n 2 matrix 
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Pab the Ky Fan norm N(pab) = Si=i °~i(PAB), which is the sum of all the singular values 
of pab, is < 1, or equivalently log(N(p~AB)) < 0, where q = min(m 2 ,n 2 ). 

Proof: Suppose pab has a decomposition of pab — J2iPiPt '® pf = Y^iPi(UfE^ ,rn 'U^)® 
(V^ E^'^V^^) with < pi < 1 satisfying J2%Pi = 1- Here is a k x I matrix which has 

A B 

entry 1 in position i,j and all other entries as zero. The U i ' are the unitary matrices which 

A B 

diagonalize p i ' . Applying the properties of Kronecker products: 

vec{XYZ) = {Z T ® X)vec(Y), (see [2S]) 

Z = X®Y Z = vec{X)vec{Y) T , (see EU) 

we have 

pf ®pf =vec{p?)vec{pf ) T 
= {Uf* ® Uf)vec(E[T' m) )((Ur ® Vf )vec{E^ n) )f 
= (Ut ® C/ i A ) U ec(^r' ro) )«ec(^' n) ) T (C/f f <8> C/f T ). 

(6) 

Since and f/j 5 are unitary, it is evident that Uf* ® and t/- 5 ^ ® Uf T are unitary 
also. Moreover, vec{E^' m ^)vec{E ( 'yl ,n ^) T — e[™ ' n ' has a unique singular value 1. Then 
Pi ® /°f nas a ^ so a unique singular value 1 due to the fact that the Ky Fan norm is unitarily 
invariant 26 . Therefore, as a convex linear combination of pf ® pf , pab should have the 
Ky Fan norm N{pab) < J2iPiN(pf <8> pf ) = = 1- Here we have used the inequality 

J2i °~i{A + B) < fi(A) + o~i(B) where A and B are fc x / matrices and q = min(fc, I) 

ESI- □ 

The cross norm criterion proposed in |23| leads to the same result as our above Theorem, 
but with the expressions in Dirac bra-ket notation. Here we obtain identical results but from 
simple matrix analysis. Moreover, our approach has special advantages in its concise and 
explicit expressions: the two factor spaces A and B just correspond to different rows and 
columns of pab under the transformation pab — y Pab ■ Also the technology for the proof is 
much simpler compared with the complicated operator algebra used in |23|. Now we only need 
to rearrange the entries of pab according to Eq. (0) , then compare the sum of the square roots 
of the eigenvalues for p^abPab^ (i- e -> the Ky Fan norm or trace norm of pab) with 1. The 
separability criterion is strong enough to be sufficient for many lower dimensional systems. 
It is straightforward to verify sufficiency for the Bell diagonal states [Ej, Werner states in 
dimension d = 2 |27j and isotropic states in arbitrary dimensions , as was done in [2.'{j . So 
far we have not yet found a close connection between this criterion^and the PPT criterion, 
but from Eq. we have the compatibility relation N(p) = N(p TA ), since p and p TA only 
differ up to a series of elementary row transformations which are naturally unitary and keep 
the Ky Fan norm invariant. For pure states, we have the following stronger conclusion: 

Proposition 1: A bipartite pure state is separable iff pab has a unique singular value 1. 

Proof: Necessity: Given a bipartite separable pure state we have pab = p A <8> p B ' ■ From 
Eq. ijBJ, Pab has a unique singular value 1. 
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Sufficiency: pab has a unique singular value o~\ = 1, so pab — <TiUiv\ = U\v\. By Eq. J3J, 
we have pab — a® (3 with vec(a) = ui and vec((3) = v\. Moreover, pab = Pab == ' > ct® ft = 
a 2 1 . Thus the eigenvalues of a and j3 should both be 1 or — 1 at the same time because 
Pab has only one eigenvalue 1. This yields p A = a and p B = (3 if a and (3 have eigenvalue 1, 
and p A — —a and p B = —ft if a and (3 have eigenvalue — 1. Therefore pab is a pure separable 
density matrix. □ 

Now we derive a dual criterion based on the Theorem: 

Corollary: For a separable to x n system pab, a permuted version pba which exchanges the 
first and second factor spaces of pab can be defined as 

Pba = S(n, m)p A BS(m, n), (7) 

where S(m, n) = Y^iiLi EjLi ■Sy™'™' ® (E^ n ' n ^) T . Then for the n 2 x m 2 realigned matrix Pba 
the Ky Fan norm N(pba) is < 1 or equivalently \og(N(pBA)) < 0. 

Proof: For a separable m x n system p^s , there exists a decomposition: pab = ^ZiViPt ® pf ■ 
Using the property Y ® X = S(n, m)(X ® Y)S(m, n) where X is an to x m matrix and F 
is an n x n one [201, we have pba = S(n,m)(J2 t Pipf <& pf)S(m,n) = Y^iPiPf ® /°f • This 
density matrix is also separable but associated with an n x to system, so we have N(pba) < 1 
according to the Theorem. □ 

This criterion is equivalent to the Theorem, due to the property Y ® X = vec(Y)vec(X) T — 
(vec(X)vec(Y) T ) T = (X®Y) T . Thus we have pa~b = (pab) T and further N(pa~b) = 
N(pba). This can also be seen from the symmetry with respect to the two subsystems 
of the cross norm given in |23j . 

We expect that log N (pab) may depict a possible measure of entanglement for the corre- 
sponding system. In fact, for the d— dimension maximally mixed state p m — Id/d 2 it is easy 
to derive pm = ^ SiJ=o NKiil ana - ^°S^{Prn) — ~ logd, while for the maximally entangled 
state p e = 2 Sfj=o l**)0'j'l we have p e = Id/d 2 and log N(p e ) = \ogd. Noticing the relations 
Pm — Pe and pm = p e , we see that they are dual and log AT (jo) is distributed symmetrically 
about where the separable pure states are located. Hence any bipartite quantum state p 
can be depicted by such an approximate measure log N(p) representing its "distance" from 
p m and p e . Furthermore we have the following property for log N(j5): 

Proposition 2: Applying a local unitary transformation leaves log AT (jo) invariant, i.e. 

log N(pYb)= log N(pa~b), (8) 

where p AB = (U ® V)pab{U^ <8> V') and U,V are unitary operators acting on the A,B sub- 
systems, respectively. 

Proof: From Eq. (J3J, we have a decomposition of pab = S;=i a i®Pi where q = min(TO 2 , n 2 ). 
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Setting ji = (U <g> V)(ai <g> ® V*), we have 

= ([/*«)C/)(wec(a 4 )wec(/3 i ))(^ t «)^ T )- (9) 

Summing all the components 7^, we have yO^ B = (U* ® U)pab{V^ ® ^ T )- Since t/ and V are 
unitary, it is evident that U* ®U and U T are unitary, hence log N{p' AB ) = log N(pab) 
since AT (jo) is a unitary invariant norm. □ 

In addition, like any unitary invariant norm, N(p) is convex |26j . but log N(pab) is not convex 
due to the concave logarithmic function. Therefore, this rough "measure" is not exact because 
the criterion is not a necessary and sufficient one, but it depicts the entanglement to some 
degree in a subtle way independent of the PPT criterion. 

3. Application of the criterion to bound entangled states 

In this section we give some typical examples to display the power of our separability 
criterion to distinguish the bound entangled states from separable states. 

Example 1: 3x3 bound entangled states constructed from unextendible product bases (UPB) 
In Bennett et al introduced a 3 x 3 inseparable bound entangled state from the 
following bases: 

hfc>> = 7fl°>(l°>- I 1 ))' IV'i) = -i=(|o> — |i»|2>, 

life) = ^IW)- |2», |Va) = ^(|l)-|2))|0), 

M = i(|0) + |l) + |2»(|0> + |l) + |2)), 

from which the density matrix can be expressed as 

1 4 

P= l (Id-Y,\^)(^\)- ( 10 ) 

i=0 

Direct computation gives log N(p) = 0.121, which shows that there is slight entanglement in 
this state. 

When the UPB have the following construction 

\i>j) = Wj) ® \v 2 j mod 5), i = 0, . . . ,4 



where the vectors Vj are Vj — TV (cos ^p, sin ^Hi, h), with j — 0, ... ,4 and h = | yl-j- \/5 and 
N = 2/\/h + v% then the PPT state of Eq. (TJJ gives logiV(p) = 0.134, which identifies 
this bound entangled state. 



Example 2: Horodecki 3x3 bound entangled state 
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Horodecki gives a very interesting weakly inseparable state in [23 which cannot be detected 
by the PPT criterion. The density matrix p is real and symmetric: 
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(11) 



where < a < 1. To compare our result with the computation for the entanglement of 
formation in |29|. we consider a mixture of the Horodecki state and the maximally mixed 
state p p = pp + (1 — p)Id/9, and show the function 

/ = max(0,logiV( / 5 p )) (12) 

in Figd 



a 




Fig. 1. Estimate of degree of entanglement for a Horodecki 3x3 bound entangled state. When 
a = 0.236 we have maximal entanglement of the state for 0.9955 < p < 1. 

We see that logiV(p^) > when p = 1 and < a < 1, which completely identifies these 
bound entangled states. Furthermore, when a — 0.236, log N(p p ) has a maximum of 0.0044 
and follows a similar trend to that of the degree of entanglement for p p (with a maximum at 
a = 0.225) given in It can be seen that we obtain an upper bound p = 0.9955 for p p 

which still has entanglement when a — 0.236. 

To test the criterion, we programmed a routine to analyze various PPT entangled states 
described in the literature. On a 600MHz desktop computer we performed a systematic search 
by checking 10 5 randomly chosen examples of the seven parameter family of PPT entangled 
states in within half an hour. The Theorem could detect entanglement in about 21% of 
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these bound entangled states satisfying p = p TA and gave a maximum value for log N(p) of 
0.22, which is a convincing demonstration of the strength of the criterion. 



Remark 1: 

One might expect the Theorem to be equivalent to the PPT criterion for a 
2x2 system. Unfortunately, this is not the case, and certain 2x2 entangled 
states cannot be recognized by our criterion. We take the following example from 
[7| which is a separable state only when p = | or a = or a = 1: 




Fig. 2. The function / = max(0, log N(p)) for the quantum state of Eg. 1131 




where a,b > and \a\ 2 + \b\ 2 = 1. We plot the value of / = max(0, log AT (jo)) in 
Fig. [21 with respect to a and p. We can see that, apart from when p = i or a = 
or 1, there are still some regions where / = 0. Thus these 2x2 entangled states 
cannot be completely detected by our criterion. 



Remark 2: 

It is interesting to make a comparison between log N(p) or N(Jj) — 1 with the 
entanglement of formation Ef |SJ. By direct comparison, for the density matrix of 
Eq-El log N {ft) or N(p) — 1 is less than the entanglement of formation calculated 
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in ■ As for the state of a 2 x 2 system, we can conveniently calculate the entan- 
glement of formation by Wootter's formula |2|. By straightforward computation, 
we have 

1. For 2 x 2 entangled Werner states (2ZJ: log AT (jo) >Ef = AT (jo) - 1; 

2. For the density matrix of Eq. EH N(p) — 1 < Ef but there is no definite 
relation between log Y(jo) and Ef, 

3. For 3 x 3 entangled isotropic states [10]: AT (jo) - 1 > E f and logAT(jo) > Ef. 

Thus we find that there is no ordered relationship between Ef and AT (jo) — 1 or 
logATGo). 



4. Summary and discussion 

Summarizing, we have presented a matrix realignment method to assess the separability 
of the density matrix for a bipartite quantum system in arbitrary dimensions. The criterion 
provides a necessary condition for separability and is quite easy to compute. It shows re- 
markable ability to recognize most known bound entangled states in the literature where the 
PPT test fails. Moreover, it gives a rough estimate for the degree of entanglement. We also 
show that this criterion is not equivalent to PPT even for 2x2 quantum states. Moreover, 
there is no definite relationship between our estimate and the entanglement of formation. 
However, in combination with the PPT criterion, we can significantly expand our ability to 
distinguish directly the entanglement and separability of any quantum state in arbitrary di- 
mensions. Our method has recently been developed further by a linear contraction approach 
that gives a generic separability criterion . It is expected that this method will find more 
applications in the study of multipartite quantum systems and other problems in quantum 
information theory and quantum computation. We leave the explicit physical meaning of the 
criterion as an open question that awaits further study. 
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